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Abstract We examined the dietary diversity and food intake of Phrynocephalus frontalis, compared the difference of 
insect diversity in the natural habitats with different lizard densities, and discussed the potential role of this lizard in the 
desert ecosystem. The results show that: (1) arthropodans of the orders Coleoptera, Hymenoptera and Hemiptera were 
major dietary components of P. frontalis; (2) coleoptera larvae always formed the predominant component of lizard 
diets; (3) dietary diversities of P. frontalis were not significantly different between summer and autumn or between 
the two sexes; (4) the similarity in trophic niches between seasons was 0.756, whereas the similarity in trophic niches 
between sexes was 0.994; (5) stomach content weight of lizards varied significantly among different seasons, but there 
was no significant difference in stomach content weight between sexes; (6) insect diversity differed significantly among 
the groups of the habitat with different degrees of lizard density, and the habitat with moderate lizards density had the 
highest insect diversity. We infer that P. frontalis prey mainly on insects and change their diet and food intake with 
season; males and females consumed similar preys in types and weights. As an important predator, P. frontalis could 


affect the insect community in the arid ecosystem of Hunshandak Desert on the Mongolian Plateau. 
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1. Introduction 


Many lizards live in arid ecosystems where productivities 
are low due to low precipitation. Thus how lizards select 
their diets and maintain food intake are critical issues 
for these lizards. When the climatic or biological factors 
restrict food intake, some lizards (e.g. Sauromalus 
obesus) change their diets, adjust the time budgets, 
and apparently defer reproductive activities to survive 
the severe environment conditions (Nagy, 1973). As 
in other animals, lizard diets also relate to intrinsic 
factors such as hunger impetus and morphological traits 
(Verwaijen et al., 2002). 
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Lizards consume large numbers of insects each year 
and assumed to play an important role in insect control 
(Knowlton, 1938; Pacala and Roughgarden, 1984). The 
latest studies have focused on the relationship between 
insect diversity and predation pressure from reptiles and 
found that lizards influence insect community structure 
(Murakami and Hirao, 2010; Cozzens, 2011; Garda et al., 
2012; Borkhataria et al., 2012). Similarly, losses of reptiles 
due to climatic changes in the Mediterranean reduce 
their potential effects on insect control (Araújo et al., 
2006; Civantos et al., 2012). 

Phrynocephalus frontalis (a common name is steppe 
toad-headed lizard) is a common sand-dwelling lizard in 
Central Asian deserts (Zhao, 1997; Munkhbaatar et al., 
2006). Previous reports showed that P. frontalis is 
a dominant predator in some sandy lands or deserts 
(Wang and Fu, 2004; Zhao, 2001; Liu et al., 2008; Li 
et al., 2011). Knowlton (1938) suggested that lizard 
impacts on insect populations depended upon diet, insect 
consumption rate and the abundance and distribution 
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of lizards. Lian et al. (2012) described the abundance, 
density and distribution pattern of P. frontalis in 
Hunshandak Desert. However, there was no quantitative 
study of the ecological trophic role of P. frontalis in the 
desert ecosystem. There, these lizards tended to aggregate 
in suitable habitat patches, and density reached to 0.174 
individuals/m’ (Lian, 2011; Lian et al., 2012). It is 
necessary to understand the potential effects of P frontalis 
on insect communities in Hunshandak Desert. 

In this study, we investigate the diet and food intake 
of P. frontalis in summer and autumn. Then, we analyze 
the influences of season and sex that influence dietary 
diversity of P. frontalis. Finally, we compare differences 
in insect diversity among habitats having different lizard 
densities and discuss the potential trophic role of P. 
frontalis in desert habitat. 


2. Materials and Methods 


2.1 Ethics statement In this study, we adhered to the 
‘Guidelines for the use of animals in research’ published 
in Animal Behaviour 1991 and the Wild Animals 
Protection Law of the People’s Republic of China. 
All animals in this study were cared for under animal 
research protocol IOZ-2006 approved by the Animal Care 
Committee of the Institute of Zoology, Chinese Academy 
of Sciences. 


2.2 Study area This study was conducted in the sand 
dunes around Sangendalai (115°57' N, 42°40' E) in 
Xilingole, Inner Mongolia, China. Sangendalai is in the 
southern region of Hunshandak Desert, with an average 
elevation of 1813 m above sea level. In this region, 
mean annual temperature is 1.7 °C, with the maximum 
temperature of 35 °C in July and the minimum of —33 
°C in January. The growing season lasts for 189 d; mean 
annual precipitation is about 300 mm, with most rain 
falling in July and August; and mean annual evaporation 
is 1936.2 mm, which is five times larger than precipitation 
(Ding et al., 2005; Nie and Zheng, 2005). The study area 
is dominated by sandy soil and also has meadows and 
saline soil in the lowland region, which is sparsely covered 
by Caragana microphylia, Artemisia halodendron, Ferula 
bungeana, Bassia dasyphylla, Agropyron cristatum 
and Cleistogenes squarrosa (Liu and Guo, 2003; 
Peng et al., 2006). 


2.3 Species identification and food intake We carried 
out observations in June (summer) and September 
(autumn) of 2011. We totally collected 56 P. frontalis 
adults (27 in summer, 29 in autumn; 31 of males, 25 of 
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females) to form our samples. The existence of hemi- 
penis was used to identify the male P. frontalis. 

During the morning (9:00-11:00 h) of each sampling 
day, we captured lizards by hand, intraperitoneally 
injected 75% ethanol solution in each individual promptly, 
and then stored lizard samples in 75% ethanol solution 
for laboratory analysis. We weighed and dissected each 
lizard, and identified the food items in their stomach 
contents to degree of Order (Tatner, 1983). We used 
stomach content weight to represent the food intake of the 
lizard. We also used Shannon-Wiener index to quantify 
dietary diversity of the lizards (Shannon and Weaver, 
1949; Fialho et al., 2000). 


H'=-> „log, P 


Where: H'= the value of the Shannon-Wiener diversity 


h . 
" species, s = the 


index, P, = the proportion of the i 
number of species in the community. 
We used the niche overlap index of Levins (Pianka, 
1973; Maia et al., 2011) to assess the degree of similarity 
of the trophic niche of P. frontalis between summer 


and autumn: 


2 2 
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Where O,, is Pianka’s measure of niche overlap index 
between two seasons j and k; P, and P, are the numerical 
proportion of prey category i in the diet in seasons j and 
k, respectively. 


2.4 Lizard density, insect diversity and plant coverage 
We chose 1520 m x 200 m transects in our study area. By 
using mark-recapture method (Kacoliris et al., 2009), we 
investigated lizard density of each transect and classified 
those 15 transects into three categories: (1) Low density, 
the habitat with low lizard density of less than 0.01 
individuals/m’; (2) Moderate density, the habitat with 
lizard density of between 0.01 and 0.1 individuals/m’; and 
(3) High density, the habitat with lizard density of more 
than 0.10 individuals/m’. 

To investigate insect diversity in the habitats, we 
arranged 20 paper cups with baits at 10 m intervals 
along the length of transect. The baits were combined 
with sugar, vinegar, 75% ethanol solution, water in the 
proportion of 1:2:1:20 in weight. We also used insect 
net to collect insects above each trap. All collected 
insects were kept in sealed plastic bag and identified in 
laboratory. We used Shannon-Wiener index to quantify 
insect diversity (Shannon and Weaver, 1949). 

In each transect, we systematically arranged ten 1 
m x 1 m wood quadrats at 20 m intervals along the 
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length of transect. Digital cameras have been used for 
remote sensing over short distances and are considered 
to be a reliable method for measuring plant coverage 
(White et al., 2000; Chi et al., 2007). By this method, 
we used Photoshop 6.0 (©1990-2002, Adobe Systems 
Incorporated) to measure plant coverage of those chosen 
1 m’ wood quadrats. The plant coverage was calculated 
with the formulas below. 


c= 100% 
A 


Where: C = plant coverage of a quadrat, G = pixels of 
green in a quadrat, A= pixels of all colors in a quadrat. 


2.5 Statistical analysis All statistical analyses were 
conducted with SPSS 13.0 (SPSS, Inc., Chicago, IL, 
U.S.A.). Kolmogorov-Smirnoy tests showed that the 
distribution of variables of dietary diversity, stomach 
content weight, insect diversity and plant coverage 
all fit a normal distribution (P >0.05). Thus, we used 
the Univariate GLM (general linear model) to test the 
difference of each variable among different experimental 
groups, as well as the interaction effects of those factors. 
If necessary, we used multiple comparison of Post Hoc 
test of Tukey to calculate the significance among pairs of 
variables. The Bonferroni correction was used for these 
multiple comparisons. Probabilities lower than 0.05 were 
considered significant for all tests. 


3. Results 


3.1 Dietary components in summer and autumn The 
lizard preyed on Insecta and Arachnida (Table 1). Plant 
remains were only found in the stomachs of lizards 
sampled in autumn. 


In summer, the first two predominant components of 
the lizard prey were the Coleoptera (with the percentage 
of 32.6% of the total weight consumed and the rate of 
occurrence of 96.3% in stomach, respectively) and the 
Hymenoptera (with the percentage of 56.39% and the 
rate of occurrence of 81.48% in stomach, respectively) 
(Table 1). 

In autumn, the predominant prey items were 
Hymenoptera (50.83% of prey items, present in 93.10% 
of stomachs), followed by Hemiptera (Table 1). Arachnida 
animals were also another important prey in autumn 
(Table 1). 

Coleoptera larvae accounted for 5.73% of prey items 
in summer and 4.79% in autumn; occurring in 37.04% of 
stomachs in summer and 44.83% in autumn (Table 1). 


3.2 Dietary diversity and stomach content weight 
Dietary diversity indices of P. frontalis were not 
significantly different between summer and autumn (F 
= 0.082, df= 1, P = 0.776; Table 2) or between the two 
sexes (F = 0.750, df= 1, P = 0.397; Table 2). The value of 
trophic niche overlap (O,,) between summer and autumn 
was 0.756, whereas the value of O, between two sexes 
was 0.994. 

Stomach content weights showed significant difference 
between summer and autumn (F = 10.358, df= 1, P = 
0.002; Table 2). The lizards preyed more in autumn than 
in summer (Table 2). There was no significant difference 
in stomach content weight between two sexes (F = 2.082, 
df= 1, P= 0.165; Table 2). 


3.3 Insect diversity, lizard density and plant coverage 
Insect diversity indices differed significantly among three 
habitats with different lizard densities (F = 4.023, df= 
2, P = 0.022; Figure 1) Post Hoc tests Tukey showed 


Table 1 Percentage and rate of occurrence of prey items in the diet of Phrynocephalus frontalis from Hunshandak Desert in Inner Mongolia 


of China, during the summer and the autumn. 


Taxon Percentage of prey items (%) Rate of occurrence (%) 
Class Order Summer Autumn Summer Autumn 
(n = 27) (n= 29 ) (n=27) (n=29) 
Insecta Coleoptera 32.16 3.54 96.3 37.93 
Hymenoptera 56.39 50.83 81.48 93.1 
Hemiptera 0.44 34.58 ae 75.86 
Lepidoptera 0.85 0.42 7.41 6.92 
Diptera 0.88 0.21 7.36 3.52 
Orthoptera 0 0.48 0 6.88 
Coleoptera (larvae) 5.73 4.79 37.04 44.83 
Lepidoptera (larvae) 0 0.42 0 3.44 
Diptera (larvae) 0 0.24 0 3.37 
Hymenoptera (larvae) 0 0.38 0 3.45 
Arachnida Araneae/Acarina 3.52 2.08 25.93 20.69 
Plants 0 0.63 0 10.34 
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that insect diversity of the habitat with moderate lizards 
density was significantly higher than that of other two 
habitats (P = 0.014 and P = 0.041, respectively; Figure 
1) There was no significant effect of plant coverage on 
insect diversity index (F = 0.860, df = 4, P = 0.492; 
Figure 2) The interaction effects of lizard density and 
plant coverage were also not significant (F = 1.596, df= 
8, P= 0.140). 


4. Discussion 


Phrynocephalus frontalis mainly preyed on Coleoptera, 
Hymenoptera and Hemiptera. The dietary diversity indices 
of P. frontalis showed no significant difference between 
summer and autumn, but the main prey items changed 
from Coleoptera in summer to Hymenoptera in autumn. 


0.8 H 
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Insect larvae were also the most important prey items of 
lizards in growing seasons (Sales et al., 2012; Taylor et al., 
2012). Our data showed that beside of imago, insect 
larvae were being increasingly preyed on by P. frontalis in 
autumn. Optimal foraging models predict that the shortage 
of more generalized diets lead predators to consume more 
“sub-optimal” prey than in periods of greater resource 
abundance (Maia et al., 2011). Thus, P. frontalis could be 
an opportunistic insectivore lizard which changed their 
prey according to the seasonally changing environment 
of Mongolian Plateau. Similar findings were reported in 
brown anoles (Anolis sagrei) that their diet varied from 
place to place and from season to season (Spiller and 
Schoener, 1990; Spiller and Schoener, 1997; Huang et al., 
2008; Norval et al., 2010). Moreover, it is suggested 
that historically shifts of dietary lead to lizard diversity 


0.4 H 


Diversity index of insects 


Oe High density 


Moderate density 
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Type of lizard density 


Figure 1 Insect diversity index in habitats having different lizard density. For the X axis, low density: 0.00-0.01 lizards/m*; moderate 


density: 0.01—0.1 lizards/m’; high density: > 0.10 lizards/m’. 
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Figure 2 Insect diversity index of the habitat with different plant coverage. For the X axis, plant coverage is: the lowest: 0-10%, lower: 
10%-20%, moderate: 20%-30%, higher: 30%—40%, and the highest: >40%, respectively. 
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Table 2 Dietary diversity and stomach content weight of Phrynocephalus frontalis and their influencing factors. * indicates the significant 


difference between different groups. 


Factors Index of dietary diversity Niche overlap index (Oj) Stomach content weight (g) 

Season Summer (n = 27) 1.377 + 0.097 0.756 0.078 + 0.009 
Autumn (n = 29) 1.373 + 0.118 0.092 + 0.009 
Significance 0.776 0.002* 

Sex Male (n = 31) 1.025 + 0.214 0.994 0.109 + 0.013 
Female (n = 25) 1.238 0.131 0.144 + 0.020 
Significance 0.397 0.165 


observed today (Vitt and Pianka, 2005). That means the 
relationships between dietary taxa and phylogeny could 
result in various dietary of Agamidae lizard in modern 
time (Vitt and Pianka, 2004; Vitt and Pianka, 2005). 

The food intake and energy demands in sand lizards 
related to the environmental changes and endocrinological 
action (Telemeco et al., 2010). For example, appetite- 
related hormone secretions of Qinghai toad-headed 
lizards (Phrynocephalus vlangalii) respond to cold 
stress slightly to mobilize energy and live their vivid 
life in Qing-Tibetan Plateau (Li et al. 2011). Base on 
our study, stomach content weights of P. frontalis varied 
significantly among different seasons. They preyed more 
on insects in autumn (non-mating season) than in summer 
(mating season). Rose (1982) found that there was no 
overt relationship between stomach content weight and 
reproductive season in the iguanid lizard Anolis acutus. 
Thus, although they belong to the same superfamily of 
Iguania, A. acutus and P. frontalis showed different prey 
consumptions in different environments they lived. 

Generally speaking, animals living under severe 
conditions would be expected to reduce their sexual 
differences in foraging patterns and food habits (Hoffman, 
1983; Caravello and Cameron, 1987). The indifference 
of food intake of male and female lizards represented the 
same food demand of both sexes. In the same way, the 
types and sizes of prey consumed by male and female 
whiptail lizards were similar (Sales et al., 2012). 

Additionally, there was a role of the lizards to their 
habitat where they could affect animal diversity by 
predation. We found that the insect diversity indices 
differed significantly among the groups of habitats with 
different lizard densities. The insect diversity of the 
quadrat with moderate density of lizards was significantly 
higher than that with high or low density of the lizards. 
Murakami and Hirao (2010) found that the presence of 
predatory lizards strongly affected species richness of 
native insects. Biologists suggested that lizards can be 
used to control invasive species for protecting native 
diversity and resisting biological invasion (Wanger et al., 


2011). Lian (2011) also reported the possibility that man 
use the P. frontalis to control pests in desert grassland. 
We assumed that P. frontalis could affect the insect 
community in the arid ecosystem of Hunshandak Desert 
on the Mongolian Plateau. 
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